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ABSTRACT 



We have conducted J, H, and Ks imaging observations for the Coalsack 
Globule 2 with the SIRIUS infrared camera on the IRSF 1.4 m telescope at SAAO, 
and determined the color excess ratio, Ej_h / Eh^Kq- The ratio is determined in 
the same photometric system as our previous study for the p Oph and Cha clouds 
O ! without any color transformation; this enables us to directly compare the near- 

^ \ infrared extinction laws among these regions. The current ratio Ej^h / Eh-Ks = 

<N ; 1.91 ± 0.01 for the extinction range 0.5 < Ej^h < 1.8 is significantly larger 

\C) \ than the ratios for the p Oph and Cha clouds [Ej^h / Eh~Ks = 1-60-1.69). This 

ratio corresponds to a large negative index a = 2.34 ± 0.01 when the wavelength 
dependence of extinction is approximated by a power law Ax oc A~", which 
might indicate little growth of dust grains, or larger abundance of dielectric non- 



O , 

^ , absorbing components such as silicates, or both in this cloud. We also confirm 

^ I that the color excess ratio for the Coalsack Globule 2 has a trend of increasing 

with decreasing optical depth, which is the same trend as the p Oph and Cha 
clouds have. 
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Introduction 



Precise determination of the wavelength dependence of extinction is necessary to un- 
derstand the properties of interstellar dust because differences in the extinction law are 
attributed to a variation in the dust grains. Differences in extinction laws among various 
lines of sight throughout the optical and UV properties of the spectrum are widely accepted 
(e.g., Fitzpatrick 1999). The existe nce of a corresponding variation at somewhat longer 
wavelengths is a controversial issue (jWhittetl Il988l : iMathid Il990l : iKenyon et al.l Il998l ) . For 
near-infrared wavelengths, recent studies (e.g., Messineo et al. 2005, Nishiyama et al. 2006) 
have found that the extinction law seems to change from one line of sight to another. These 
results indicate that the so-called universality (e.g., Mathis 1990) does not necessarily hold 
for infrared wavelengths. 

The color excess ratio, Ej^h / Eh-Ksi is one of the simplest parameters available to 
express the near-infrared extinction law, which can be determined only from J, H, and 
Ks photometric observations. In the observations, however, different photometric systems 
and transformation between them can cause difficult problems in comparing the extinction 
laws in different regions. Naoi et al. (2006) have shown that the photometry transformation 
between different photometric systems does not reproduce the real color excess ratio probably 
because the transformation is mainly determined from the photometry of intrinsically red 
stars, not reddened stars. Hence, Naoi et al. (2006) have conducted imaging observations 
for the p Oph and Cha clouds with the same photometric system, and have not found any 
significant differenc es in the extinction law, Ej^h I Eh-Kq , between the clouds in contrast 



to previous works (IKenyon et al.l Il998l : iGomez fc KenyonI l200ll ). Naoi et al. (2006) have 
suspected that the apparent difference of the color excess ratio, Ej_h / Eh-Ks^ was caused 
by the transformation between different photometric systems, at least for these two dark 
clouds. 

Racca et al. (2002) observed the Coalsack and determined the color excess ratio 
Ej_h/Eh-k = 2.08 ± 0.03, which is significantly larger than the ratios for other clouds 
(the mean value of 1.61 ± 0.04, Whittet 1988). They suggested possible changes in grain 
chemistry with extinction or star formation activity. Though the color excess ratios for the 
p Oph (1.57 ± 0.03, Kenyon et al. 1998), Cha (1.80 ± 0.03, Gomez & Kenyon 2001), and 
Coalsack (Racca et al. 2002) clouds were all expressed in the same CIT system, the pho- 
tometric system used in the actual observations are different, and the color transformation 
between them made the comparison very difficult. 

Hence we have conducted J, H, and Ks imaging observations for the Coalsack Globule 
2 (about 3 mag deeper than that of Racca et al. 2002) with the same photometric system 
(SIRIUS/IRSF) as in Naoi et al. (2006) for the p Oph and Cha clouds, and determined 
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the color excess ratio Ej^h / Eh-Ks- We discu ss the extinctioii law of the Coalsack Globule 
2 in addition to the p Oph and Cha clouds (INaoi et al.l l2006l ). and also discuss the grain 
properties in comparison with extinction model calculations. 



2. Observations 

We acquired J (1.25 /im), H (1.63 yum), and Ks (2.14 /im) images of the Coalsack 
Globule 2 in March-April 2004 with the Simultaneous Infrared Imager for Unbiased Sur- 
vey (SIRIUS) on the Infrared Survey Facility (IRSF) 1.4 m telescope at the South African 
Astronomical Observatory (SAAO) in Sutherland, the Republic of South Africa. Dichroic 
mirrors enabled simultaneous observations in the three bands (see Nagashima et al. 1999 
and Nagayama et al. 2003 for the details on the instrument). The camera is equipped with 
three 1024x1024 pixel HgCdTe (HAWAII) arrays. The pixel scale of the array was 0'.'453 
pixel"^, giving a field of view of 7f7x7'7. 

We acquired 20 s exposure sets for 3 times. Ten dithered frames were obtained for one 
set of exposures. Total integration time is 10 minutes for each position. Relative positional 
offsets of dithering were 20"or 25". Typical seeing conditions were 1'.'3, 1'.'2, and I'.'l (FWHM) 
in the J , and Ks bands, respectively. The observations were made at air masses between 
1.0 and 1.5. Dark and twilight flat frames were obtained at the beginning and end of each 
night. We observed standard stars from the faint near-infrared standard-star catalog of 
Persson et al. (1998) for photometric calibration on the same night. 

The target region encompasses heavily extinguish ed part of the Globule 2 in Coalsac k 



determined from optical extinction and radio CO maps (IRacca et al.ll2002l : iKato et al.lll999l ). 



We also observed a reference region in order to estimate the intrinsic co lor, which lies ver y 
close to the target region and falls in regions of relatively low extinction (IRacca et al.ll2002l ). 
The observed areas are summarized in Table [H 

We used NOAO's Imaging Reduction and Analysis Facility (IRAF)E| software package to 
reduce the data. We applied the standard procedures of near-infrared array image reduction, 
including dark current subtraction, sky subtraction, and fiat fielding. Each image, following 
subtraction of the average dark frame, was divided by the normalized fiat-field image. Then 
the thermal emission pattern, the fringe pattern due to OH emission, and the reset anomaly 



^IRAF is distributed by the National Optical Astronomy Observatory (NOAO), which is operated by the 
Association of Universities for Research in Astronomy (AURA), Inc., under cooperative agreement with the 
National Science Foundation. 
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slope pattern of the HAWAII arrays were corrected for each frame, with subtraction of 
a median sky frame. Identification and photometry of point sources in the all frames were 
performed with using the DAOPHOT packages in IRAF. The lOcx limiting magnitudes (< 0.1 
mag) for point sources were J ~ 19.0, H ~ 18.3, and Ks ~ 17.0 mag. We used these point 
sources in analysis, but the point sources brighter than J < 12.0, H < 11.0, and Ks < 10.0 
(less than 1% of the total) were not included in our samples because they were saturated. 



3. Results and determining of the color excess ratio Ej^h / Eh-Ks 

Figure [T] shows J — H vs. H — Ks color-color diagrams of point sources for the Globule 
2 and the reference region in Coalsack. The photometry data described in this paper is in 
the SIRIUS/IRSF photometric system because we do not use any transformation equations 
to other photom etric systems for photometry in order to avoid possible sources of errors 



( iNaoi et al.ll2006l ). In Tabled], the numbers of identified point sources are also shown. 



In order to determine the color excess ratio, Ej^h / Eh-Ks^ fo^' the Coalsack Globule 2 
we employ the powerful method developed by Kenyon et al. (1998). The method determines 



the r atio accurately even when a number of PMS stars are present in the sample (INaoi et al. 



20061 ). We assume that the stellar population behind the cloud is identical to the stellar 
population in the reference region. The point source distribution in a color-color diagram 
is transformed to the density distribution by a kernel density estimator with the kernel as 
triweight. Figure [2] shows J — H vs. H — Ks density distributions drawn from Figured] for 
the Globule 2 and the reference region in Coalsack. Contours are drawn in a logarithmic 
scale (the same as for all figures below). These color distributions use a kernel density 
estimator with a smoothing parameter h = 0.2 to derive the density, which is roughly 
twice the 1 a error of our photometry at the survey limits. The sum of multiplications of 
the density distributions between the target and reference regions provides the reddening 
probability distribution, which is defined as the probability of measuring a pair of color 
excesses, Ej_h and E^-Ks- Figure [3] shows Ej^h vs. Eh-Ks reddening probability contours 
for the Coalsack Globule 2. We derive the slopes of the reddening vector, namely the color 
excess ratio Ej^h / Eh-Ksi by tracing the ridge of contours on the reddening probability 
distribution which is divided into the circular annuli centered on the origin. The slope of 
the reddening locus Ej_h/Eh-Ks^s determined by a least-squares fit of a straight line with 
data having errors in both x and y coordinates (see Naoi et al. 2006 for more details). The 
analysis yields Ej_h/Eh-Ks = 1-91 ± 0.01 (for 0.50 < Ej_h < 1-8) in the SIRIUS/IRSF 
photometric system, and the result is summarized in Table [2] In Figure [S] the enlargement 
with plots of the ridgeline on contour and the best fit lines for the different extinguished 
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parts are also shown for clearness. 



4. Extinction law in the Coalsack Globule 2 

Table [2] shows the compilation of the color excess ratios, Ej h / Eh-Ko. for the Coalsack 



Globule 2 and for the p Oph and Cha clouds (INaoi et al.ll2006l ). The color excess ratio for 
the Coalsack Globule 2 (1.91 ± 0.01) is significantly larger than those for the mean values of 
the p Oph (1.65 ± 0.03) and Cha (1.67 ± 0.01) clouds. Moreover, the value in the Coalsack 
Globule 2 is higher than those in any sub-regions of the p Oph and Cha clouds. These results 
indicate that a real variation in the near-infrared extinction law occurs from cloud to cloud 
(see also, Nishiyama et al. 2006 for the case of the observations toward the Galactic Center). 
The result in the Coalsack Globule 2 supports the relatively large ratio indicated by Racca 
et al. (2002). The shape of the extinction curve is frequently characterized by a law of the 
form A(^X) A~", where a is the spectral index for near-infrared extinction (e.g., Draine 
2003). A good fit to the curve for J, H, and Ks ba nds is obtained with a = 2.34 ± 0.01 



by using the mean effective wavelength of the filters (iNishiyama et al.ll2006l ): this index is 
very large compared to the values in other regions (e.g., Whittet 1988 shows a mean value 
of a = 1.70 ±0.08). 

Such a large color excess ratio, Ej^h/ Eh~Ks^ ^ large spectral index, a, in the 
Coalsack are realized when the upper cutoff of the grain size distribution of radius is kept 
around 0.2 pm, and also if we assume greater dielectric material ratio. Figure H] shows the 
calculations of the color excess ratio and the index a as a function of the upper cutoff of 
the grain size distribution of radius, a+, and the silicate/graphite ratio in the simple MRN 
model in Naoi et al. (2006). This may indicate that dust grains in the globule do not grow in 
size so much compared with the standard grain size ~ 0.1 pm in diffuse interstellar medium 
(e.g., Mathis 1977). Also, similar model calculations indicate that non-absorbing dielectric 
materials such as silicates cause rapidly decreasing extinction towards the longer wavelength, 
and result in large color excess ratios and large spectral indices because the efficiency of the 
scattering cross section varies in inverse proportion to the fourth power of the wavelength 
(Rayleigh scattering, e.g., Bohren & Huffman 1983). 

We also confirm a systematic shift in the major axis of the probability contour in Figure 
El The plots of the ridgeline on contour and the best fit lines for the different extinguished 
parts are also shown in the enlargement. The color excess ratio has a trend of decreasing 
from Ej^h/Eh-Ks = 2.02-2.08 for 0.5 < Ej^h < 1-0 (solid line) to Ej_h/Eh-Ks = 1-87- 
1.91 for 1.3 < Ej^H < 1-8 (dashed line). The trend is consistent with those in Naoi et 
al. (2006) for the p Oph and Cha clouds. Kenyon et al. (1998) suggests that the color 
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excess ratio in the p Oph cloud decreases from Ej^h / Eh-k = 1-65 for 0.3 < Ej^h < 0.4 to 
Ej_h/ Eh-k = 1.56-1.57 for 1.5 < Ej_h < 1-8. Kenyon et al. (1998) also indicates a large 
color excess ratio Ej_h/ Eh-k ~ 2 for the much weaker extinction range at Ej_h < 0.2 
in the cloud. They suspect that this change might be statistically insignificant. In our 
case, however, the larger number of reference sources allows us to determine the ratio more 
precisely. In order to verify the trend, the color excess ratios are determined for the two 
groups of the sample areas which are defined by the distance from the center of the Globule 
2 in Kato et al. (1999); (a) the "inside" sample (radius < 250 arcsec), and (b) the "outside" 
sample (radius > 250 arcsec). The results are complied in Table [2l The color excess ratio for 
the "inside" sample Ej_h/Eh-Ks{^'^) = 2.01 ± 0.03 is smaller than that for the "outside" 
sample Ej^h / En-Ksiout) = 2.50±0.05 in the same optical depth range of 0.7 < Ej^h < 1-2. 
The trend of the systematic shift in the major axis of the probability contour in Figure [3] 
is thus verified. As suggested by Naoi et al. (2006) for the p Oph and Cha clouds, this 
indicates grain growth with optical depth, toward the inside of the dark cloud, also in the 
Coalsack Globule 2. Therefore, although the grain size in the Coalsack Globule 2 seems to 
be small compared with other dark clouds, a slight change in size might be taking place from 
its periphery to the central region. 



5. Conclusions 

We have conducted J, H, and Ks band imaging observations for the Coalsack Globule 2 
with SIRIUS/IRSF 1.4m at SAAO and determined the color excess ratio, Ej_H/EH_Kg, 
by using the method in Kenyon et al. (1998) and Naoi et al. (2006). 

The color excess ratio for the Coalsack Globule 2 has a s ignificantly larger value than 



those for the p Oph and Cha clouds (INaoi et al.l l2006l ). Because of employing the 



same photometric system, our studies are simple and straightforward in comparison of 
near-infrared extinction laws. 

3. The extinction law for the Coalsack Globule 2 shows a trend that the color excess ratio 
decreases with increasing optical depth. The trend indicated grain growth inside of the 
cloud, consistent with the results in Naoi et al. (2006) for the p Oph and Cha clouds. 

4. The large color excess ratio for the Coalsack Globule 2 could be attributed to small 
grain growth and grain composition differences. 



We are grateful to the SIRIUS/IRSF and SAAO staff, Yoshifusa Ita, Noriyuki Mat- 
sunaga, Yoshikazu Nakada, Toshihiko Tanabe, and Cynthia Strydom for their help in support 
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of observations. T.N. are also grateful to Koji Kawabata for the calculation of extinction, 
and University of Tokyo, NAOJ, and ISAS/JAXA staffs for support and assistance in this 
study. We wish to thank the anonymous referee and editors for careful comments. This 
work is supported by grants-in-aid from the Japanese Ministry of Education, Science, and 
Culture (Nos. 12309010, 16340061, 16077204, and 15340061). 
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Fig. 1. — J — H vs. H — Ks color-color diagrams of point sources for the Globule 2 and the 
reference region in Coalsack. 
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Fig. 2. — J — H vs. H — Ks density distributions drawn from Fig{T] for the Globule 2 and 
the reference region in Coalsack. 
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Fig. 3. — Ej_H vs. Eh-Ks reddening probability contour for the Coalsack Globule 2 and 
the enlargement. The plots of the ridgehne on contour (cross marks) and the best fit hues 
Ej^H = 2.05£;i/_/f5 - 0.06 (foi" 0-5 < Ej^h < 1-0; solid line) and Ej^h = 1.89Eh-Ks + 0M 
(for 1.3 < Ej^H < 1-8; dashed line) are also shown in the enlargement. The bold parts of 
the hues correspond to the fit ranges of the ridgeline. 
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Fig. 4. — Model calculations of the color excess ratio, Ej^h / Eh-Ks (solid curves) and 
the spectral index, a (dashed curves) as a function of the upper cutoff of the grain size 
distribution of radius, a+, for the different silicate/graphite ratios in MRN (53/47; bold 
curves) and 90 % (90/10; fine curves). The large color excess ratio and spectral index in the 
Coalsack can be produced by increasing the abundance ratio of the silicate grain which is 
the representative material of non-absorbing grains. Spectral index a is fitted for the mean 
effective wavelengths o f the J, if, and Ks bands in the SIRIUS/IRSF photometric system 
dNishivama et allbooeh . 
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Table 1: Center coordinates (J2000.0) of the target and reference regions and the number of 
identified point sources. 



Region 


R.A., Dec. 


Number 


Globule 2 


12 31 24, -63 45 11 


7236 




(io:o X lo'.o) 




Reference 


12 33 20, -64 24 39 


15679 




(12:0 X 12f0) 





Table 2: Determined color excess ratios, Ej^h/ Eh-Ksi for the Coalsack Globule 2 with the 
p Oph and Cha regions^ (SIRIUS/IRSF photometric system). 



Region 


Ej-h/Eh-Ks 


Range 


(Ej-h) 


Coalsack 


Globule 2 


1.91 ±0.01 


0.50- 


-1.80 




inside*^ 


2.01 ±0.03 


0.70- 


-1.20 




outside^ 


2.50 ±0.05 


0.70- 


-1.20 


p Oph 


L1688 


1.60 ±0.01 


0.50- 


-2.35 




L1689 


1.66 ±0.01 


0.50- 


-1.70 




L1712 


1.68 ±0.01 


0.50- 


-1.50 


Cha 


I 


1.69 ±0.01 


0.50 


-2.00 




II 


1.66 ±0.01 


0.50- 


-1.80 




III 


1.67 ±0.02 


0.50- 


-1.30 



"The color excess ratios Ej-h/Eh-Ks^ot the p Oph and Cha clouds are from Naoi et al. (2006) 

''The samples of the "inside" group are selected by the distance (radius < 250 arcsec, the number of point 

sources is 3248) from the center coordinate in Kato et al. (1999). The samples of the "outside" are all the 
other point sources except for the "inside" samples (the number of point sources is 3988). 



